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Time-periodic driving provides a promising route to engineer non-trivial states in quantum many- 
body systems. However, while it has been shown that the dynamics of integrable, non-interacting 
systems can synchronize with the driving into a non-trivial periodic motion, generic non-integrable 
systems are expected to heat up until they display a trivial infinite-temperature behavior. In this 
paper we show that a quasi-periodic time evolution over many periods can also emerge in weakly 
interacting systems, with a clear separation of the timescales for synchronization and the eventual 
approach of the inhnite-temperature state. This behavior is the analogue of prethermalization 
in quenched systems. The synchronized state can be described using a macroscopic number of 
approximate constants of motion. We corroborate these findings with numerical simulations for the 
driven Hubbard model. 
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I. INTRODUCTION 

Experiments with ultra-cold atomic gases in optical 
lattices and ultra-fast spectroscopy nowadays allow to 
address the dynamics of quantum many-particle systems 
out of equilibrium. A particularly important role in this 
context is played by periodically driven systems m- 
Periodic driving can stabilize novel states both in cold 
atoms and in condensed matter, including topologically 
nontrivial states [5HIQ] or complex phases such as super¬ 
conductivity mm- It can be used to engineer artificial 
gauge helds in cold atoms m and emergent many-body 
interactions such as magnetic exchange interactions in 
solids pAtilB] . or to transiently modify lattice structures 
through anharmonic couplings m- 

An important question is thus the theoretical under¬ 
standing of the long-time dynamics of periodically driven 
systems. The approach to a steady state has been investi¬ 
gated intensively for the relaxation of isolated systems af¬ 
ter a sudden perturbation, both experimentally and the¬ 
oretically [TBl - ED] . When a generic non-integrable many- 
body system is left to evolve with a time-independent 
Hamiltonian, it is believed to eventually relax to a ther¬ 
mal equilibrium state, unless it is in a many-body lo¬ 
calized phase pT1123j . If the system is integrable, on the 
other hand, the steady state is often described by a gener¬ 
alized Gibbs ensemble (GGE) [MIES], which keeps track 
of a macroscopic number of constants of motion. When 
integrability is only slightly broken, the system can dis¬ 
play dynamics on separate timescales, such that observ¬ 
ables rapidly prethermalize to a quasi-steady nonequilib¬ 
rium state which can be understood by a GGE based 
on perturbatively constructed constants of motion |26j . 
before thermalizing on much longer time scales p7fl29j . 

Integrability turns out to be a crucial factor also for 
periodically driven systems. Their dynamics can syn¬ 
chronize with the driving |30j and display a non-trivial 


periodic time evolution at long times. A way to under¬ 
stand this is to show that the time evolution over one 
period T commutes with an infinite number of opera¬ 
tors Xx, which are thus conserved at stroboscopic times 
(i.e. integer multiples of the period). Having a fixed ex¬ 
pectation value of all Xx at stroboscopic times, one can 
construct a statistical ensemble to describe the long-time 
behavior of the system (the periodic Gibbs ensemble), 
which has been analytically and numerically shown to 
give correct predictions for hard-core bosons En¬ 
in contrast to integrable systems (and many-body lo¬ 
calized states [32H37] ). it has been proposed that generic 
non-integrable systems “heat up” under the effect of driv¬ 
ing and display rather trivial infinite temperature proper¬ 
ties as soon as they settle into a periodic motion [55HI0] . 
One can formulate this statement in terms of the Eloquet 
eigenstates (the exact solutions of the Schrodinger equa¬ 
tion with a periodic evolution of all observables 0111]), 
stating that each individual Eloquet state displays infi¬ 
nite temperature properties. This conjecture relies on 
a breakdown of the perturbative expansion of Eloquet 
eigenstates at some order because of unavoidable reso¬ 
nances between transitions in the many-body spectrum 
with multiples of the driving frequency. A common ap¬ 
proach to avoid this problem is to construct effective 
Eloquet Hamiltonians from a high-frequency expansion 
[ISEl]. In this work we show that a quasi-periodic state 
can also emerge in weakly interacting systems, provided 
that linear absorption can be avoided: the stroboscopic 
time evolution is constrained by approximately conserved 
constants of motion Xx- Analogously to prethermaliza¬ 
tion in weakly interacting systems after a sudden pertur¬ 
bation the system rapidly synchronizes with the 

driving and remains periodic over a large number of peri¬ 
ods m such that 3> mT ^ J“^, where g controls the 
strength of the interaction, T is the period of the driv¬ 
ing and J is coupling of the integrable Hamiltonian, e.g.. 


2 


the bandwidth of the kinetic energy term; we set h = 1 
throughout. The quasi-periodic state can be described 
as a periodic Gibbs ensemble based on the I\, i.e., stro¬ 
boscopic prethermalization gives access to quasi-periodic 
states which are entirely different from the infinite tem¬ 
perature final states. 

This paper is organized as follows: in Sec.|TT]we intro¬ 
duce the general formalism: first, we rotate the weakly 
interacting Hamiltonian in such a way that it commutes 
with the integrals of motion of the noninteracting part 
(Sec. IIA I, then in Sec. IIB we identify the approximate 
integrals of motion and finally in Sec. |II C] study the time 
evolution of the observables. In Sec. IIIII we discuss the 
relation with the periodic Gibbs ensemble |3Ij and in 
Sec. |IV| the relation with the Floquet theory of periodi¬ 
cally driven systems. In Sec.|V]we specialize the results of 
Sec.|ll]to the Hubbard model, first presenting our analyt¬ 
ical results (Sec. V A | and then comparing them to the 


numerical findings (Sec. V B) obtained with dynamical 
mean-field theory. In Sec. VI we draw our conclusions. 


II. GENERAL FORMALISM 


In the following we consider an integrable, noninter¬ 
acting system perturbed by a weak periodic interaction. 
The general Hamiltonian is given by 

H{t)=Ho+gH,^t{t), (I) 

where the integrable part 

Ho=Y,exIx. (2) 

A 

can be written as a sum of constants of motion (e.g., 
momentum occupations for independent particles on a 
lattice), and the small parameter g controls the strength 
of the interaction iJint(t) which is periodic with period 
T and frequency H = 27:jT. To study the time evolu¬ 
tion at stroboscopic times tm = rnT (to integer), we ex¬ 
tend the approach of Refs. [26l [28] to periodically driven 
systems, and determine a time-periodic unitary trans¬ 
formation R{t) such that the Hamiltonian in the 

rotated frame commutes with the constants of motion 
I\ at any time up to corrections or order 0[g^) [42]. If 
\^{t)) = R{t)\ip{t)) is the transformed wave function, the 
Hamiltonian Fleff which dictates the evolution in the ro¬ 
tated frame via idt\'4’{t)) = Ffeff(i)IV’(^)) is given by 

H^s{t) = R{t)H{t)R{t)'< - iR{t)R{t)'f . (3) 


Since R{t) is assumed to be unitary, we make the ansatz 
R{t) = with an anti-hermitian operator S{t). 

In the next sections, we first analytically find the trans¬ 


formation S(t) (Sec. HA), then identify the approximate 
integrals of motion (Sec. HB[) and finally in Sec. H C 


we 


study the time evolution of the expectation values of ob¬ 
servables and the dependence of their long-time behavior 
on the frequency H. 


A. The transformation S{t) 


We now show in detail how to rotate the Hamiltonian 
Eq. 0 with a transformation R{t) = such that the 
Hamiltonian 0 is (i) periodic and (ii) diagonal in the 
operators that diagonalize Hq. To implement condition 
(i), we first expand to second order in g and then write 
in Fourier series both the effective Hamiltonian 


Hesit) = H!:l\t) + gHi^\t) + g^H^^it) 

E —inflt r ^(0) _L ,^2 ^(2) 


( 4 ) 


and the anti-hermitian operator 

n 


( 5 ) 


with HeS,n = Hlff and S'„ = -Sln. Combining 
Eqs. 0 and 0 we find: 

Hesit) = i?o + <? (^i?int(t) + [S^^Ht),Ho]+if^S^^Ht)^ 
+ 5" (^l[S^^\t),Ho] + [A«( 0 ,ffint(t)] + 
^-[SW^t),[S^^\t),Ho]] + lj^S^^Ht)+ 
-^(5(i)(t)5«(t)-^W(t)5(i)(t))^ +0(g3) . 

( 6 ) 


To ensure condition (ii), we require that 

[HSnM=0, (7) 

for any Fourier component, perturbative order and con¬ 
stant of motion, labeled by n, X and A respectively. As 
in Ref. 1^ we employ the basis Tx\a.) = aAlc^)- We as¬ 
sume that the energies ex are incommensurate, so that 
the eigenenergies of i7o, he. Ea = are non¬ 

degenerate. For an extensive lattice model this can be 
achieved, e.g., by using sufficiently irregular boundaries. 

After some lengthy but otherwise straightforward alge¬ 
bra, we can find Hf.s(t) and S(t) by repeatedly applying 
Eq. Q to each perturbative order in Eq. 0 , so as to 
reduce the Hamiltonian to the diagonal form 

HMt) = + l«)^diag,c.(0(a| + 0{g^) . (8) 


In order g'^ we have 


(0) _ J i7o 


if n = 0 
otherwise , 


( 9 ) 


so that \a){a\E''^>^, with E'^^l = E, 


(0) 


AO) 
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To first order in g the Fourier components of S{t) read: 

( 10 ) 


( (/Slgint.nla) jr U- R 

0 otherwise 


The first-order perturbative correction to H^s is: 


where 




E^ui = (a|i?int.„|a) . 


a 


( 11 ) 


( 12 ) 


At order the Fourier components of are found 
to be: 




0 ( 1 ) TT I Zjy^ 

Jp , -flint,n-p -I- -eidiag,n- 


( 1 ) 


|a) 


— Ea — TlVL 


(13) 

if a ^ /3 and, as previously, we choose the diago¬ 
nal elements to be zero. In Eq. (13) we have defined: 


-ffdiig.n = \oi)Eli'}a{a\. Finally, the second order term 
of the effective Hamiltonian reads: 




(14) 


with: 




(a|gint,p|/?)(^|gint,n-p|a) 
Ea- Ep- pfl 


Pita p 

- (°^l-^int,n-p|/3) (/3|-gint,p|«) 

Ep — Ed — pH 


(15) 


B. Approximate integrals of motion 

Under a general unitary transformation, the time prop¬ 
agator U(t,0) = Te“*lo ) is transformed into 


[/(t, 0 ) tj{t,0)e^^°\ 


(16) 


with U (f, 0) = Te“® lo ^ Because S{t) is periodic, 

the time evolution at stroboscopic times is thus unitar- 
ily equivalent to the time evolution with the diagonal 
Hamiltonian ([s]), 

tmO{g^). This implies that the quantities 


Xx = 


(17) 


are approximately conserved under the evolution over 
multiple periods T, i.e., {ix{tm)) = (Ia( 0 )) + tmO{g^). 
For the example of a weakly interacting Hubbard model 
studied below, the original constants of motion are mo¬ 
mentum occupations rife of independent particles, while 
the constants of motion of the stroboscopic time evolu¬ 
tion correspond to quasiparticle modes. 


C. Expectation value of observables 

We examine the synchronization of these modes in 
terms of the time evolution 


(A)t = (^(0)|C7l(t,0)AC/(t,0)|V>(0)) 


(18) 


of an observable A which is a function of the original con¬ 
stants of motion Xx (having in mind, e.g., a measurement 
of momentum occupations rife or higher-order momentum 
correlation functions nfen^/), assuming that the system 
is in an eigenstate |'0(O)) = |a) of EIq before the driving 
is switched on. Inserting Eq. (16) into (18), expanding 
the operators and in powers of g, and using 

the fact that [A,U{t,0)] = tO{g^) (because A commutes 
with all Xx), we obtain 

(A}t = -2Re(a|S'(0)A[S'(0) - 5(t)]|a) -b tO(g^), (19) 

with S(t) = U\t,0)S{t)U{t,0) . For stroboscopic times, 
with S{tm) = *S'( 0 ) determined by Eq. ( 10 ), one finds the 


final result for the perturbative time evolution 


(A),, = E 


^2 / ojt„ 




)ynp{yj) 


".p_r 


(w — nO^jiuj — pQ) 


+ t^O{g^), ( 20 ) 


where ?/np(w) denotes the spectral density 

ynp{to) ={cx\H{dt^—nA6(^uj Hq + Ed)H[Yit,—p\A) . ( 21 ) 


The integral in Eq. (20) gives an accurate description of 
{A)t^ for times tm ^ g~^, where relative corrections 
traO{g^) are small. Note that for finite m the term 
sin(tma;/2) regularizes the singularities at nuj. There¬ 
fore the amount of contributing spectral weight is due to 
the location of H inside or outside the band, as discussed 
below. For 5 —>■ 0 there is thus a large time window 
g~^ ^ tm ^ T in which the dynamics is governed by the 
long time asymptotics of the integral. To analyze this, 
we distinguish two different behaviors depending on the 
frequency H: 

(i) Fermi golden rule regime: If there is nonzero spec¬ 
tral density y„„(nH) > 0 at an even pole l/(w — nH)^, 
the stroboscopic evolution for m 1 develops a linear 
asymptotics {A)t^ ~ Z)/ 3 /a(/ 3 |^|/ 3 )ra->/ 3 , where 

is the Fermi golden rule excitation rate 

F„^^ = 2^^ |(/?|77i„t,„|a)P<5(nH - Ep + E^). (22) 


To see this fact one can consider the contribution to the 


integral ( 20 ) from a small interval \uj — nH| < e around 
the pole, in which ynn(w) can be approximated by a con¬ 
stant ynninU). With a substitution x = tm{oj — nfl), the 
remaining integral is tm dx sir?{x/ 2) /x'^ ~ tmT^l2,. 
From a similar consideration for n A P ocie can obtain 
the subleading terms. 

(a) Stroboscopic prethermalization: Assuming that the 
perturbation involves only a limited number of Fourier 
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components, such as for a harmonic perturbation with 
Hint,n = 0 for |n| > 1, then the spectral density ynm{i^) 
is restricted to a finite band [—VF], depending on 
the type of excitation, the bandwidth of the noninter¬ 
acting single-particle spectrum, and phase space restric¬ 
tions. If all poles uj = nVl lie outside this band, the limit 
TO —)■ oo integral of (201 is simply obtained by replacing 
sin^(tmw/2) by its average 1/2, which corresponds to the 


first term in Eq. (191 


(A)p,e = -2Re(a|^(0)Fl^(0)|a) 


=vi: 


duj 


ynp(w) 


n,p 


(uj — nf})(uj — pfi) 


(23) 


In this case the system synchronizes for ^ T (and 
tm ^ 9~^) into a periodic evolution with values {A)t^ = 
(A)pre, before further heating takes place on longer 
timescales. This is the analogue of prethermalization in 
a quenched system. 


In principle, one can now use standard first-order pertur¬ 
bation theory to construct perturbative Floquet states 

\i’a,n) = H-, where the zeroth order 

is given by the unperturbed eigenstates = ^m,o|ct) 

(^Fq = Ea). The perturbative expansion does not con¬ 
verge to the true Floquet eigenstate if there are reso¬ 
nances Ea — Ejj = nfl in the many-body spectrum, but 
low orders nevertheless can exist: in particular, the first 
order is given by and it is well-defined 

outside the Fermi golden rule regime. This shows that 
the prethermalized state Eq. (23) is related to the per¬ 
turbative Floquet state by 


(4l)p.e = 2(V-).,L|dI|^W ). 


^ F^Oil 


(26) 


Here the appearance of a factor of two is reminiscent to a 
similar relation between the prethermalized and ground 
state expectation values in the quench case. 


III. STATISTICAL DESCRIPTION OF THE 
PRETHERMALIZED STATE 


V. APPLICATION TO THE HUBBARD MODEL 
IN INFINITE DIMENSIONS 


The condition j/„„(nr2) = 0 for the absence of linear 
absorption is equivalent to the absence of resonances in 


Eq. (10). Outside the Fermi golden rule regime, the con¬ 
stants of motion 0 are thus well-defined and one can 
ask whether the prethermalized state can be described by 
a Gibbs ensemble pQ = I^g EH: where the 

Lagrange multipliers px are determined by the constraint 
from the initial state, (Ia)o = 


(A)pre = tr[pgA]. 


(24) 


Using Eqs. (23) and ( [l7| ), the proof for this statement 
only relies on the time-independent matrix S'(O) being 
antihermitian and appearing only to order g in tr , 
and thus proceeds analogously to the argument showing 
that prethermalized states for a sudden quench can be 
described by a GGE [26] . 


IV. RELATION TO THE FLOQUET PICTURE 


We now explain how the prethermalized state Eq. (23) 


can be related to the Floquet spectrum of the Hamil¬ 
tonian. According to the Floquet theorem, the ex¬ 
act solution of the Schrodinger equation with a time- 
periodic Hamiltonian (IH is given in the form \'4’F,a{t)) = 
g-zFy,at where^^/Q(t)) is periodic in time. If a 

system is in a Floquet state, the time evolution of ob¬ 
servables is periodic. By expanding \%l^a(t)) in a Fourier 
series |t/'a(i)) = the Floquet quasi¬ 

energy spectrum can be obtained by diagonalizing the 
time-independent block-matrix, 

(L/F,a T toQ E{]i')\'lpa,m) ~ 9 ^^, Ejnt^l |'0a,m+/) ■ (25) 


A. Analytical results 

In order to illustrate the general results above, we now 
choose as specihc example the Hubbard model 

H{t) = -J cl^Cja + U{t)Y^{n^^-^){n^i-l), (27) 

{ij)a i 

with nearest neighbor hopping J and periodically mod¬ 
ulated interaction 


U{t) = 17(1 — cos(Qt)). 


(28) 


With these choices, the first and the second term of 


Eq. (27) represent the integrable, noninteracting part i7o 


and the periodic weakly interacting perturbation with 
g = U, respectively. Energy and time are measured in 
units of J and J~^, respectively. The constants of motion 
of Hq are momentum occupation numbers fik^ = cl.^Cka- 
To allow for a comparison of the analytical results de¬ 
rived above and a numerical solution, we consider the 
model in the limit of infinite spatial dimensions with a 
semi-elliptic density of states p{e) = a/ 4 — e^/(27r) at half¬ 
filling (density n = I). In this limit, the dynamics can 
be computed using nonequilibrium dynamical-mean-field 
theory [33| , and iterative perturbation theory [151SS] as 
impurity solver (see Sec. VB). 


To investigate the prethermalization dynamics we use 
the momentum occupations as observables, A{t) = 
n(efc, t) — n(efe, 0), where n(efc,0) is the initial occupation 
of the single-particle state with energy e^. For the har¬ 
monic driving in Eq. (28) we have Hmt.n = h„i7int with 
Hint = - i)(»^4 - i)i and ho = 1, h±i = -4. 

We now proceed as in Refs. [15] and noting that 
in those derivations also initial free thermal states are 
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allowed by virtue of the finite-temperature version of 
Wick’s theorem. The time-dependent occupation of a 
state with single-particle energy e at time tm = fnT is 
given by [ 221 SS] : 

^pert(Cfc, tm) — 0) dtt ^ra) i (^9) 

where we choose the initial distribution 7 T,(efc, 0 ) = 
(cj[.^c/ccr)t=o to be thermal, and 


F(e,tm)=^ fduj ^ 


{ui — nfl){u} — pfl) 


(30) 


where we have dropped the fc-dependence since momen¬ 
tum conservation can be omitted in the limit of infinite 
dimensions (i.e. one has Jfc(w) = We find: 


Je(a;) = / deide2de3 p{ei)p{e2)p{e3)[n{e3)n{e)n{ei)n{e2) 


- n{ei)n{e 2 )n{e 3 )n{e)]S{ei -|- £2 - £3 - w - e), (31) 

where h(e) = 1 — n(e) (which equals n{—e) in the case 
of particle-hole symmetry, which we consider here). The 
function Je(a;) (Eq. ( |M| )) has already been obtained for 
the investigation of the sudden quench [ 22 ] (which is 
contained in our results by setting h±i = 0). The 
connection with the spectral density (211 is given by 

2/e,np(^) — 


From Eq. (311, one can read off the phase space condi¬ 


tion for the Fermi golden rule: at zero temperature, n(e) 
= 0 (—e) and p(e) = 0 for |e| > 2 , hence linear absorption 
(Jfc(±r 2 ) 7 ^ 0 ) should occur for for |efc| < 11 < 6 -|- |efc|. 
More details on the phase-space argument leading to ei¬ 
ther the Fermi-golden rule regime or the prethermaliza- 
tion plateau can be found in the Appendix, where useful 
expressions for the numerical evaluation of Eq. (31) are 
also presented. 


B. Numerical results 


In Fig. we show the single-particle occupation n(efc) 
at stroboscopic times for a specific value of e for H = 3.93 
and 12 = 10.47, which lie in the Fermi golden-rule regime 
and in the prethermalization regime, respectively. We 
find that the perturbative predictions from Eqs. (20) 


and (31) capture well the initial slope of the occupation 


in the linear absorption regime, as well as the prethermal¬ 
ization plateau predicted by Eq. (23) for 12 = 10.47. For 


later times the numerical results approach the infinite- 
temperature value n^_>o(e) = 0.5. As expected, the 
agreement between the DMFT results and the perturba¬ 
tive predictions improves with decreasing U, where the 
prethermalization plateau extends to longer times. In 
the inset of Fig. [^ we show the time evolution of the oc¬ 
cupation n(e) at U = 0.8. At 2 ~ 2-4 the quasi-periodic 



FIG. 1. Momentum occupation n(efc) for energy = —0.4 
(upper panels) and double occupation (lower panels) at stro¬ 
boscopic times 27rm/12 for U = 0.1, 0.3, 0.5, 0.8,1.2 from top 
to bottom curves in each panel. Panels (a) and (c): 12 = 3.93, 
panels (b) and (d): 12 = 10.47. Symbols with dashed lines: 
DMFT data, continuous lines: perturbative predictions from 
Eqs. (201 and (311. Inset of panel (b): same as main panel 
for U = 0.8, but showing complete time evolution (i.e. also 
non stroboscopic times). The arrows in panel (b) show the 
predictions of Eq. (|23[). 


prethermalization regime begins where n(e) is constant at 
stroboscopic times. The double occupation (lower pan¬ 
els of Fig. also shows a prethermalization plateau at 
high frequency, while it evolves towards its infinite tem¬ 
perature value for f2 in the Fermi golden-rule regime. In 
Fig. i we plot n{e,tm) as a function of e after a given 
number of periods (m = 11). Panel (a) corresponds to a 
frequency such that every value of e gives rise to linear 
terms, which are on the contrary absent for Vl = 10.47, 
see panel (c). Panel (b) refers to an intermediate case 
(T = 1.0, 12 = 27r), where only the boundary values of 
e (i.e. e > —2 and e 2 ) give linear contributions and 
thus at tm ^ T differ from the DMFT data. Finally, 
panel (d) shows the absorption of energy, measured by 
the slope am(e, 12) = n(e, m27r/12) — n(e, (m — l)27r/12), 
which becomes small for 12 > 6 — e, as predicted by the 
perturbative calculation (shown with a dashed line for 
m = 11). We point out that the regime of validity of 
the DMFT calculation with iterative perturbation theory 
does not allow to explore small values of the frequency 
(12 < 1 ) where the other boundary (e = — 12 ) lies. 


VI. CONCLUSIONS 


In conclusion, we discussed the analogue of prethermal¬ 
ization in periodically driven systems. A weakly interact¬ 
ing system can synchronize into a quasi-steady state with 
nontrivial properties, before reaching the infinite temper- 
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for the specific case of the driven Hubbard interaction by 
rewriting Eq. (30) and applying a phase-space argument. 


As a first step, we express Je(w) in terms of 


OO 

R{s) = j de n{e) p{e) , (Al) 


using a Fourier representation of the delta function: 




R{sf. 

(A2) 


FIG. 2. (a), (b), (c): single particle occupations for a driv¬ 
ing term with U = 0.5, computed with DMFT (red triangles, 
with initial inverse temperature j3 = 20) and from perturba¬ 
tion theory (blue dashed lines, initially in the ground state) at 
an intermediate time, tn = IIT, for Q. — 27r/r = 3.93, 6.28, 
and 10.47, respectively; the initial state is prepared with U 
= 0. The green continuous line shows the free initial ther¬ 
mal state at (3 = 20. The UPT prediction in panel (c) neg¬ 
ligibly differs from the prethermalization plateau predicted 
by Eq. ( |23| |. (d): Excitation over one period, measured as 
am(£, H) (see main text) for t/ = 0.1 using DMFT. The black 
dashed line corresponds to e = —fl + 6, which together with 
e = — H delimits the Fermi-golden-rule regime. 


ature state generic for the long-time behavior of driven 
non-integrable systems. This stroboscopic prethermal¬ 
ization is a consequence of the existence of a macroscopic 
set of operators which are almost conserved by the time 
evolution over one period. Stroboscopic prethermaliza¬ 
tion thus provides a way to engineer quantum states with 
a nontrivial effective dynamics, alternative to the a high 
frequency expansion. These states reflect the properties 
of perturbative Floquet states, which can be very differ¬ 
ent in nature from the exact Floquet states. 
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We also note that for an initial zero-temperature state 
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Fn,p\^i im) — 


(n — 

K.n(e, tm) = (e, im, nn ), 


(A3) 


where we defined 


OO 

FW(e,t„,F;)= J 


dui 




UJ 


N 


J,iuj + E). (A4) 


Consider first the case of zero (or sufficiently low) tem¬ 
perature of the initial state and |e| < |r2| < 3D + \e\. Then 
a term linear in tm contributes to F(e, tm), namely (E = 
|nH|, N = 1,2, X = {uj - E)tm) 




Trf 

6n2^ME) (A5) 
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Appendix A: Phase-space argument for the 
Hubbard model in infinite-dimensions 


As discussed in Sec. m the single-particle occupa¬ 
tions Eq. (29) at long times (i.e. tm T) display two 
regimes, namely the Fermi-golden rule absorption regime 
and the stroboscopic prethermalization regime, depend¬ 
ing on the value of and e. Here we discuss these regimes 


This corresponds to the Fermi golden rule regime with 
a linear-in-time growth of n(e,tm)- On the other hand, 
if H is outside the indicated interval, the denominators 
are never zero (for zero temperature) and a stroboscopic 
prethermalization plateau is attained. 

In all cases we can rewrite the integrals more compactly 
by using the identities 


sin^(a;t/2) 

LU 


OO 

f s\T?(ujt/2) 


cos(a;s) 


0 


- y du sin(a;u), 

0 

^(t-s)0(f-s). 


(A6) 

(A7) 
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and taking the symmetries of the w and s integrals into 


account. We obtain 


tm 

Risf X (A8) 

0 

tm 


- / dslm 




1 

2 


J dsRe 


R{sf X 


(A9) 


0 

- h(e)e-*(^+^)") . 

These expressions are suitable for numerical evaluation; 
they can be further simplified for the zero-temperature 
case. 
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